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ABSTRACT 


^  A  program  to  investigate  the  characteristics  of  integrated  optical 
components,  coupling  between  optical  waveguides  and  other  optical 
components  and  waveguide  fabrication  technology  was  initiated  in 
February  20,  1 98 l  for  Naval  Avionics  Center  under  contract  N00I63  8l  M  1798 
The  work  performed  under  this  contract  is  summarized  in  this  report.  The 
Navy  technical  monitor  was  Dr.  K.  J.  Jones. 

A  prism  coupling  measurement  system  was  modified  and  improved. 

Numerical  techniques  for  computing  the  refractive  index  and  the  thickness 
of  the  optical  films  from  the  experimental  data  are  discussed.  Techniques 
based  on  optical  fibers  and  TV  cameras  for  measuring  the  waveguide  attenua¬ 
tion,  and  techniques  related  to  end-fire  coupling  and  near  infrared 
radiation  are  considered.  Methods  for  fabricating  solution-deposited 
waveguides  and  waveguides  on  glass,  lithium  niobate  and  lithium  tantalate 


substrates  are  also  reviewed. 
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1.0  INTRODUCTION 


Since  the  useful  bandwidth  of  a  coherent  source  is  proportional  to 

its  oscillating  frequency,  and  the  frequency  of  a  laser  radiating  in  the 

\k 

visible  or  near  infrared  spectra  is  of  the  order  of  10  Hz,  its  useful 
bandwidth  is  tremendous.  Many  optical  devices  and  systems  have  been 
conceived  to  take  advantage  of  the  tremendous  bandwidth  of  lasers.  However, 
the  weight,  size  and  alignment  problems  of  bulk  optical  devices  and 
their  susceptibility  to  mechanical  vibrations,  air  turbulence  and  tempera¬ 
ture  fluctuation  prevent  these  devices  from  becoming  practical.  Recently 
interest  has  turned  to  integrated  optical  devices  and  systems,  since  i 

many  of  the  problems  associated  with  bulk  optical  devices  are  absent  in 
the  integrated  or  guided-wave  devices. 

In  February  1981,  an  integrated  optical  device  program  was  initiated 
at  Technology  Associates,  Inc.  for  Naval  Avionics  Center  (NAC) .  The 
objectives  were  to  provide  the  technical  support  for  the  enhancement  of 
NAC's  in-house  integrated  optics  (10)  activities  in  upgrading  its  experi¬ 
mental  facilities  and  evaluation  capabilities,  and  to  assist  NAC  personnel 
in  screening,  evaluating  and  assessing  appropriate  10  material,  components 
and  technologies  for  specific  systems  of  current  and  future  interst  to 
NAC. 

Optical  waveguides  are  the  basic  building  blocks  of  any  10  system. 
Optical  components  can  be  built  from,  and  optical  subsystems  can  be 
connected  by  the  optical  waveguides.  The  potential  exists  that  all 
optical  components  and  associated  electronics  may  be  fabricated  on  a 


common  substrate. 
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Waves  can  be  coupled  Into  or  out  of  the  waveguides  via  prism,  end* 
fire  or  grating  coupling.  So  far  as  experimental  investigations  are  concern¬ 
ed,  prism  couplers  are  probably  the  most  versatile  tool  for  coupling 
radiation  into  or  out  of  waveguides  and  for  characterizing  guided  modes. 

On  the  other  hand,  end-fire  coupling  is  the  most  reliable.  The  prism 
(Figure  l)  and  end-fire  (Figure  2)  coupling  schemes  will  be  discussed 
in  Sections  2  and  3  respectively.  Grating  coupling,  not  attempted 
in  this  work,  will  not  be  discussed  in  this  report.  The  effective  index 
of  refraction  and  the  attenuation  constant  are  two  basic  and  important 
characteristics  of  optical  waveguides.  Techniques  to  estimate  these  wave¬ 
guide  parameters  are  discussed  in  Sections  4  and  5.  Experiments  involv¬ 
ing  infrared  or  near  infrared  radiation  are  discussed  in  Section  6. 

Methods  of  fabricating  practical  optical  waveguides  are  presented  in 
Section  7- 


2.0  PRISM  COUPLING 

A  dielectric  slab  waveguide  consisting  of  a  substrate  region  with 

an  index  of  refraction  ns,  a  film  region  of  an  index  n^,  and  a  cover 

region  with  an  index  nc,  is  the  simplest  possible  optical  waveguide 

structure.  Many  modes  can  be  supported  by  the  structure,  depending  upon  the 

free-space  wavelength  A,  the  thickness  h  of  the  firm,  and  the  values 

of  n  ,  nr  and  n  .  Let  N  rr  be  the  effective  index  of  refraction  of  the 
s’  f  c  eff,m 

m-th  quided  mode,  then  N  --  is  related  to  these  parameters  via  the  disper- 
J  7  eft,m 

si  on  relation  [ I ,  2] : 


(I) 


-3- 


2irh  /  2 
1T/nf 


^  rr  *  V  (n^,n  ,n  ,N  „  ) 

eff#m  m  r  s  c  eff,m 


where 


VWWf.m*  "  ""  +  ♦s(nrVM«ff,»)  +  *c(nf ’VNeff  J  (2) 


♦.(nf'n5'Maff(.)  =  arctan 


♦c(nf*nc»Neff,m)  "  arctan 


[m^>] 


(3) 


W 


and  p  *  0  for  TE  modes  and  p  *  1  for  TM  modes. 


2.1  Basic  principle  of  prism  coupling: 

If  a  prism  with  an  index  n^  and  an  angle  e  is  placed  next  to  the  film, 
(see  inserts  of  Figures  3  and  4)  and  an  optical  beam  is  incident  upon  the 
boundary  between  the  film  and  the  prism  at  an  angle  such  that  the  projection 
of  the  phase  velocity  of  the  incoming  beam  in  the  prism  matches  exactly  with 
c/Neff >  e  being  the  speed  of  light  in  vacuum,  of  the  m-th  guided  mode,  then 
the  m-th  guided  mode  will  be  excited.  This  is  the  basic  principle  of  prism 
coupling.  The  angle  of  incidence  a  or  a",  n^e  and  are  related  by  the 

relations  below  [l,  2]: 


2  2  1/2 

*  sinefn^  -  sin  a)  +  cosc  sina 


M  .  /  2  .2  ^ 1/2  .  , 

m  sine^n^  -  sin  a  )  -  cose  sina 


(5) 


(6) 


-4- 

The  angles  a,  a'  and  e  are  defined  in  the  inserts  to  Figures  3  and  4. 

2.2  Design  considerations  for  prism  couplers: 

The  key  component  of  a  prism  coupler  is,  of  course,  the  prism. 

The  following  points  should  be  taken  into  account  in  selecting  prisms: 

(1)  The  prism  should  be  as  small  as  practical. 

(2)  The  prism  material  should  be  as  hard  as  possible  to  minimize 
the  effects  of  scratching  and  chipping  [2], 

(3)  The  index  of  the  prism  should  be  somewhat  larger,  not  merely 
slightjy  larger,  than  the  index  of  the  film,  e.g.,  n^  >  n^.  +  0.1  or  n^  >  n^.  +  0.2 

(4)  To  avoid  spurious  response  due  to  repeated  reflections  from  the 
prism  surfaces,  one  should  avoid  prisms  with  angles  of  30°,  45°,  60°,  90° 
or  other  rational  fractions  of  360°  [2]. 

(5)  If  a  half-prism,  as  depicted  in  Figures  l,  3  and  4,  is  used, 
the  angle  at  the  corner  where  the  guided  modes  are  launched  should  be 
sharp  and  well  defined. 

(6)  The  prism  should  have  a  flat  top,  rather  than  a  sharp  apex, 
so  that  pressure  can  be  applied  uniformly  to  the  prism. 

A  possible  design  of  prism  is  shown  in  Figure  5.  For  a  flint  glass 
prism  (with  n^  =*  1.75)  of  the  shape  shown  in  Figure  5,  the  dependence  of 
Neff  as  runct'ons  a  or  calculated  from  (5)  or  (6),  are  shown  in 
Figure  3.  If  the  prism  is  made  of  an  uniaxial  crystal  like  rutile  (TiC^) 
with  its  optic  axis  parallel  to  the  edge  of  91.0°,  i.e.,  perpendicular  to 
the  figure  shown,  the  extraordinary  index  ne  *  2.87  should  be  used  for  np. 


'4,  '  * 
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Figure  4  depicts  the  relationship  between  and  a  or  a"  for  such  a 

rutile  prism.  Theoretically,  coupling  to  substrates  with  values  ranging 

from  1.4179  to  2.7036  can  be  achieved  with  the  rutile  prism  and  these  values 
correspond  to  the  coupling  angles  of  a"  =*  90.°  and  a  »  90.°  respectively. 

In  reality,  the  useful  range  is  limited  to  1.7295  (a'  *  40. °)  to  2.6112 
(a  =  50.°)  region. 

2.3  Operations  of  prism  couplers: 

In  addition  to  a  finite  number  of  guided  modes,  a  waveguide  structure 
can  support  a  continuum  of  radiation  modes  in  the  cover  (air)  region  above 
the  film  and  the  substrate  region  below  the  film.  These  fields  are  referred 
to  as  the  air  and  substrate  modes  respectively.  The  air  mode  has  an 
oscillatory  field  distributions  in  all  regions  and  has  an  effective  index 
of  refraction  smaller  than  n^.  When  Neff  is  increased  beyond  n^,  there 

is  the  substrate  mode  which  has  an  exponentially  decaying  field  distribution 
in  the  cover  region  and  oscillatory  field  distributions  in  the  film  and 
substrate  regions.  For  the  guided  modes,  n^  <  <  n^,  and  the  field  decays 

exponentially  in  the  air  as  well  as  the  substrate  region. 

As  a  specific  example,  consider  a  waveguide  consisting  of  a  Corning 
7059  glass  film  (n^.  *  1.555)  of  thickness  h  *  1.0  pm  and  a  substrate  of  pyrex 
glass  (ns  =  1.45),  of  thickness  of  0.1  cm.  So  far  as  the  guided  modes  are 
concerned,  the  pyrex  glass  can  be  considered  as  infinitely  thick.  Two  TE 
modes  may  be  supported  by  the  Corning  7059/pyrex  structure  with  =  1.5356, 

and  1.4780.  Suppose  a  glass  prism  with  an  index  of  n^  =  1.75  as  shown  in 
Figure  5  is  used  in  the  experiment.  For  a"  >  27.°,  the  air  mode  is  excited. 
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A I  though  the  air  mode  is  reflected  from  the  substrate-air  boundary, 
the  reflection  coefficient  is  small  (|p|  *  0.034).  If  the  laser  beam 

is  not  blocked  by  any  objects,  it  would  emerge  behind  the  pyrex  glass. 

As  the  incident  angle  is  increased  with  a'  <  27.°,  or  a  <  10. 4°,  the 
substrate  mode  is  excited.  The  optical  beam  is  totally  reflected  at  the 
air-substrate  boundaries.  As  a  result,  the  beam  is  trapped  in  the  pyrex 
glass.  A  dot  appears  at  each  point  where  the  reflection  occurs.  Since 
the  reflection  is  very  efficient,  with  |p|  =*  I .  ,  a  series  of  dots  become 

visible.  The  distance  between  dots  is  a  function  of  the  thickness  of 
substrate  and  will  vary  with  a.  For  a  pyrex  glass  of  thickness  of  0.1  cm, 
the  distance  between  dots  on  the  same  boundary  would  vary  from  0.10  cm  to 
1.27  cm  as  the  angle  changes  from  =  25.°  to  a  =  10.0°.  As  a  approaches 
13*6°,  a  continuous  streak  would  appear,  and  this  is  the  TEj  mode  with 

=  1.4780.  When  a  is  increased  to  20.°,  TEQ  mode  with  Neff  =  1.5356 
would  be  exci ted . 

From  the  above  discussion,  three  points  helpful  to  the  proper 
interpretation  of,  and  the  successful  operation  of  the  prism  coupling 
experiments  become  obvious: 

(l)  If  possible,  there  should  be  no  reflective  object  behind  the 
substrate.  Thus  the  air  mode,  if  it  is  excited,  will  not  be  reflected. 

If  this  cannot  be  done,  then,  the  space  behind  the  substrate  should  be 
made  accessible.  By  probing  the  space  behind  the  substrate,  one  can 
differentiate  the  air  mode  reflected  by  the  object  behind  the  substrate 
from  the  substrate  mode  reflected  by  the  boundary  of  the  substrate. 
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(2)  From  the  coupling  angle  where  a  series  of  dots  begins  to  appear, 
one  can  use  (5)  or  (6)  to  estimate  the  index  of  refraction  of  the  substrate. 

(3)  As  the  coupling  angle  is  changed,  the  distance  between  dots 
increases.  The  limiting  case  should  be  very  close  to  the  highest  guided 
mode. 


3.0  END-FIRE  COUPLING 

As  mentioned  previously,  end-fire  coupling,  also  known  as  butt  coupling, 
is  one  of  three  methods  available  for  coupling  light  into  or  out  of  an 
optical  waveguides.  As  implied  by  the  name,  light  is  coupled  into  the 
waveguide  through  the  end  surface  of  the  waveguide,  as  shown  in  Figure  2. 

In  any  optical  waveguide,  energy  is  concentrated  in  a  thin  layer,  with  a 
thickness  of  the  order  of  a  few  wavelengths  at  most  near  the  waveguide 
surface .  The  region  near  the  corner  of  the  waveguide  is  critical  to  end-fire 
coupling.  Techniques  have  been  developed  specifically  to  polish  the  corner 
region  of  the  substrate  [33.  In  connection  with  this  work,  a  jig  for 
polishing  the  end  surfaces  of  optical  waveguides  has  been  designed  and  is 
to  be  used  in  conjunction  with  the  existing  polishing  facility  at  NAC. 

This  work  is  currently  in  progress. 


3.1  Transformation  of  Gaussian  beams  by  lens  or  lenses; 

To  achieve  maximum  coupling  into  the  optical  waveguide,  the  cross 
section  and  beam  divergence  of  the  incoming  optical  beam  must  match  that 
of  the  waveguide.  The  optical  beam  emitted  by  a  HeNe  laser  is  typically 
of  the  Gaussian  form.  In  particular  the  irradiance  distribution  of  a 
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TEHQq  Gaussian  beam  is  given  by  X^e  /a>  .  The  beam  radius  w  Is  defined 

-2 

as  the  radius  where  the  Irradiance  is  reduced  to  e  of  its  peak  value^Q 

at  the  beam  center  (r  *  0) .  The  beam  radius  l>q  (or  diameter)  at  the  beam 

waist,  i.e.,  the  point  with  the  smallest  beam  radius  (or  diameter),  is 

usually  specified.  In  the  Fraunhofer  zone,  the  full  angle  of  divergence 

is  2X/(ttwo).  For  example,  a  Spectra  Physics  Model  145  laser  has  a  beam 

diameter,  of  0.5  mm  and  the  beam  diverges  at  full  angle  of  1.7  mrad.  For 

these  lasers,  the  beam  waist  is  usually  located  at  the  output  mirror. 

A  simple  lens  with  a  focal  length  of  f  can  be  used  to  change  the  beam  radius. 

Let  a)  and  w  .  be  the  radii  of  the  beam  waists  of  the  incoming  and  out- 
oa  ob 

going  beams  as  shown  in  Figure  6.  woa»  anc*  their  locations  zQ  and  z^ 
with  respect  to  the  lens  are  related  by 


,  _  Z  2  7TU)  2 

_!_  «  JL  m  — i)  +  ( — 2±) 

2  2  U  f  '  V  fx  ' 

U>  .  0) 

ob  oa 


zb =  f  + 


(za  -  or 


(z  -  f)  +  (irui  /x)‘ 
a  oa 


where  X  is  the  wavelength  of  the  radiation  [4].  If  the  beam  waist  cj 

oa 

of  the  incoming  beam  is  located  at  the  front  focal  point  of  the  l<*ns,  i.e., 
z^  *  f,  then  the  beam  waist  u)ob  is  located  at  the  back  focal  point  z^  *  f 
and  has  a  value 


Ob  7TU) 


For  example,  a  lens  with  a  focal  length  of  10  cm  wi  1 1  transform  a  HeNe 


laser  beam  (A  *  0.633  ym)  with  u  of  .25  mm  to  a  beam  with  u>  .  =  80.59  ym. 

oa  od 

Since  the  cross  section  of  a  planar  waveguide  is  thin  and  wide,  it 
may  be  necessary  to  use  a  cylindrical  lens  to  transform  the  circular 
optical  beam  into  an  elliptical  shape.  The  above  equations  may  be  used 
twice,  one  for  the  major  axis  and  one  for  the  minor  axis. 

The  beam  collected  by  a  microscope  objective  or  a  lens  may  be  too 
wide  for  some  application.  Two  lenses,  with  focal  lengths  f^  and  f^  and 
spaced  at  a  distance  f^  +  may  be  used  to  change  the  beam  diameter 
by  a  factor  of  f 2/^ i >  as  shown  in  Figure  1.  By  choosing  K  the 
beam  diameter  is  reduced.  If  f^  >  fj,  on  the  other  hand,  the  beam  is 
expanded.  If  a  small  pinhole  is  inserted  at  the  common  focal  point  of 
these  two  lenses,  the  lens-pinhole  combination  also  acts  as  a  spatial 
f i 1  ter. 


3.2  Horns  for  channel  waveguides: 

Typically,  a  channel  waveguide  is  3  to  6  ym  wide  and  I  to  2  ym  thick. 
Thus  channel  waveguides  are  thin  and  narrow.  Although  the  prism  coupling 
or  the  end-fire  coupling  may  be  used  to  couple  light  into  a  channel  wave¬ 
guide,  it  is  difficult  to  align  the  optical  beam  with  respect  to  the 
channel  waveguide.  To  relieve  the  alignment  difficulty,  a  horn  structure 
may  be  used  to  transform  a  planar  waveguide  to  a  channel  waveguide  gradually. 
Linearly,  exponentially  or  parabol i cal ly  tapered  horns  have  been  proposed 
by  various  authors  [5  “  9l-  As  explained  by  Burns,  Hilton  and  Lee  [8], 
the  most  practical  horn  is  of  the  parabolic  shape.  For  parabolic  horns, 
the  width  W  of  the  horn  is  given  by 

W  -  /2kX  2  +  W 
9  o 


(10) 
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where  k  is  a  constant  and  X^  ■  X/N is  the  local  effective  guide  wave¬ 
length  of  the  channel  waveguide.  Wq  and  other  parameters  are  as  shown  in 
Figure  7.  Unless  the  channel  width,  W,  is  very  narrow  such  that  the  lowest 
order  mode  (in  the  x-direction  as  shown  in  Figure  7)  is  near  to  its  cutoff 
point,  Xg  can  be  approximated  by  XQ/n^  [8].  As  a  specific  example,  again 
consider  a  HeNe  laser  with  a  beam  diameter  of  0.5  mm  (u >  -  250  pm)  at  its 

beam  waist.  Two  lenses  with  fj  =  |6  f£  and  spaced  at  a  distance  of  fj  +  f£ 
(Figure  1)  can  be  used  to  reduce  the  beam  radius  to  15.6  pm.  The  beam 
diameter  at  the  entrance  point  is  31*2  pm.  Therefore  Wmax  =  kO  pm  is 
chosen.  Horns  have  been  designed  to  transform  W  ■  40  pm  to  channel 

iTiaX 

waveguides  of  width  5,  10  and  20  ym  respectively.  Suppose  that  the  channel 
waveguides  and  the  tapering  horns  are  to  be  built  on  photoresist  waveguides 
with  n^  *  1.6.  With  the  selection  of  k  ■  0.8,  all  parameters  of  (10)  are 
specified.  The  equations  describing  the  horns  for  5*  10  and  20  ym 
channel  waveguides  are  respectively: 

(11) 

(12) 

(13) 

For  z  *  0,  W  ■  5  ym,  10  ym  and  20  ym  respectively,  and  W  =  40  ym  when  z 
exceeds  the  values  specified  in  (II),  (12)  or  (13). 


W  *  25.0  +  0.6328  |z| ,  0  <_  z  <_  2488.9  ym 


WZ  »  100,0  +  0.6328  | z |  0  <  z  <  2369.6  ym 


W  =  400.  +  0.6328  | z | ,  0  <  z  <  1895.7  ym 


4.0  FILM  INDEX  AND  THICKNESS 


One  of  the  most  useful  features  of  prism  coupling  is  the  relationship 
between  a  (or  a**)  and  N^^.  When  a  or  a'  is  measured  experimentally,  and 
when  n^  and  e  are  known,  N^^  can  be  computed  from  (5)  and  (6).  This 
feature,  in  conjunction  with  the  dispersion  relationship  (l),  can  be  used 


to  estimate  and  h. 

For  a  waveguide,  ng ,  n^  and  X  are  usually  known  while  and  h  are 
unknown.  If  one  of  the  coupling  angles  is  measured  and  is  computed 
from  (5)  or  (6),  we  have  one  equation  ((!))  yet  two  unknowns.  There  is 


insufficient  information  to  solve  for  n^  and  h.  When  two  or  more  coupling 

angles  are  determined,  the  index  n^  and  thickness  h  of  the  thin  film  can 

be  evaluated,  provided  that  the  indices  of  the  substrate  and  the  cover, 

n  and  n  ,  are  known.  Cases  where  three  or  more  waveguide  modes  are 
s  c 

determined  have  to  be  treated  differently  from  the  cases  where  two  wave* 
guide  modes  are  determined.  In  either  case,  considerable  numerical  computa 
tlon  is  needed  to  deduce  the  final  results  from  the  experimentally  measured 


data.  A  computer  program  written  to  carry  out  the  necessary  iteration  and 
computation  is  listed  in  the  Appendix.  The  theoretical  foundation  and 


the  algorithm  of  the  computer  program  are  explained  in  this  section.  Also 
included  are  typical  results  obtained  by  the  computer  program.  The  input  and 
output  instructions  are  contained  in  the  program  MAIN.  Additional  functions 
of  the  program  MAIN  and  various  subroutines  will  be  explained  in  the  following 
sections. 


4.1  Waveguides  with  two  modes  measured: 

Suppose  that  two  coupling  angles  are  determined  and  N^^  m  with  m  *  p 
and  v  are  computed.  Then,  substituting  Neff  and  ^eff  v  'nt0  H)  anc* 
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rearranging  the  resultant  equations,  we  obtain 


2 

nf 


(N  ff  )2V  2 
eff,u  v 


<wv 


y  2  _  v  2 
U  V 


(14) 


This  is  a  nonlinear  equation  for  n^  and  a  simple  iteration  scheme  can  be 
used  to  solve  for  n^.  Specifically,  a  trial  value  for  n^  is  assumed  and 
(14)  is  used  to  compute  a  new  n^.  If  the  new  n^.  differs  from  the  old  n^. 
considerably,  the  new  n^  is  used  to  compute  a  newer  n^.  This  process  is 
repeated  until  a  satisfactory  accuracy  is  achieved.  Once  n^  is  known,  it 
is  a  simple  matter  to  solve  for  h  from  (1).  The  computation  for 
from  a  is  carried  out  at  the  beginning  of  the  program  MAIN.  The  second 
part  of  program  MAIN  is  related  to  the  iteration  for  n^  and  the  simple 
calculation  for  h. 


4.2  Waveguides  with  three  or  more  modes  measured: 

The  cases  where  three  or  more  waveguide  modes  are  measured,  assuming 
the  waveguide  has  three  or  more  modes,  are  much  more  complicated.  Suppose 
that  M  coupling  angles  are  measured,  with  M  3.  Each  measured  coupling 
angle  corresponds  to  a  value  for  m,  which  in  turn  leads  to  an  equation 

with  a  specific  value  for  m*  There  are  M  equations  (M  3)  and  yet 

there  are  only  two  unknowns.  One  cannot  solve  for  the  unknowns  from  these 
equations  unless  M-2  equations  are  redundant.  In  the  ideal  situation  where 
the  experimental  results  are  infinitely  precise  and  without  error,  M-2 
equations  will  be  redundant.  In  any  real  experiment,  however,  there  will  be 


experimental  error  or  inaccuracy;  n^  and  h  cannot  be  determined  from  the 


measured  coupling  angles.  Instead  of  "solving"  for  n*  and  h  from  N  cc  , 

t  err ,m 

the  best  estimates  for  n^  and  h  are  sought.  Reasonable  values  of  n^,  h 

are  assumed  and  used  in  conjunction  with  the  known  values  of  ng  and  nc  to 

compute  the  effective  index  of  refraction  from  the  dispersion  relation 

(1).  The  effective  index  of  refraction  so  obtained  is  labeled  as 

N  cc  (nr*h).  The  values  of  nr  and  h  are  then  searched  so  that  N  rr 
err,m  r  r  et  t  ,m 

match  with  m(n^,h).  Following  Ulrich  and  Torge  [2],  an  error  sum 

is  defined  as: 

a(nf,h)  »^Neff>.  -  Nefffi(nf,h)]2  (15) 

If  the  experimental  results  are  perfectly  accurate,  we  should  look  for  n^ 
and  h  such  that  o  =  0.  Because  of  the  presence  of  the  experimental  error 
and  inaccuracy,  this  is  not  possible.  Instead  we  look  for  the  values  of  n^ 
and  h  which  minimize  a,  and  they  are  found  by  requiring 


(16) 


9a 

3h 


■  0 


(17) 


A  gradient  method  [2]  is  used  to  solve  for  n^.  and  h  from  the  simultaneous 
equations  (16)  and  (17),  and  this  is  done  by  calling  the  subroutine  GRDT. 

The  basic  idea  of  the  gradient  method  is  quite  simple.  Suppose  a  relief 
or  a  surface  of  a  is  plotted  as  a  function  of  (n^h).  Let  a  at  a  particular 


point  (n^,h)  be  evaluated.  If  a  ■  0,  this  particular  point  corresponds  to 
the  desired  solution.  If  a  >  0,  then  a  point  with  a  smaller  value  of 
o  should  be  chosen.  If  one  is  standing  on  the  relief,  then  the  best  possible 
move  is  to  move  in  the  direction  of  steepest  descent,  which  is  precisely  the 
direction  of  the  negative  gradient.  In  GRDT,  0  for  a  given  (n^,h),  and 
its  eight  neighboring  points,  are  calculated  by  repeatedly  calling  the  sub¬ 
routine  ERSUM.  In  the  process,  it  is  necessary  to  solve  for  the  effective 
index  of  refraction  from  (1)  numerically.  A  quadratic  interpolation 
algorithm  expounded  by  Muller  is  used  for  this  purpose  [10].  To  make  use 
of  the  MULLER  subroutine^,  the  dispersion  relation  (l)  is  coded  as  the  sub¬ 
routine  DISP.  A  valley  is  reached  when  (16)  and  (17)  are  satisfied. 

4.3  Examples: 

To  test  the  accuracy  of  the  program,  two  test  runs  have  been  made. 

The  parameters  of  the  waveguide  and  prism,  the  “observed"  coupling  angles 
and  the  results  of  the  computation  are  listed  below. 

(I)  Waveguide  1: 

n  =  1.5000,  n  =  1.0000,  n  =  1.7552,  e  *  50.0°. 
s  c  p 

n^.  =  1.6500,  h  =  2.0000  pm, 
a  *  35. 95°,  32.80°,  28.80°. 

Computed  n^  *  1.6503,  Error  *  0.018%. 

Computed  h  «  2.0352  pm,  Error  *  1.76%. 


•15“ 


(2)  Waveguide  2: 

n$  -  1.4500,  nc  -  1.0000,  np  -  1.7552,  e  ■=  50.0°. 

nf  -  1.5500,  h  -  2.8000  pm, 

a  -  21.10°,  19-75°,  17.10°. 

Computed  *  1.5495,  Error  *  0.03%. 

Computed  h  -  2.7348  ym,  Error  *  2.3%. 

The  coupling  angles  are  read  and  interpolated  from  a  curve  similar  to 
Figure  3  or  4  and  are  accurate  to  ±0.1°.  Comparison  of  the  computed  values 
for  and  h  with  those  assumed  values  shows  that  the  results  are  quite 
good,  despite  of  the  inaccuracy  in  re&Mng  the  coupling  angles. 

5.0  ATTENUATION 

Although  the  attenuation  of  the  guided  beam  can  be  determined  by 
using  two  or  three  prisms  [II,  12],  the  data  so  obtained  is  not  very 
reliable.  Instead,  an  optical  fiber  may  be  used  to  probe  the  intensity 
of  the  guided  beam  and  the  attenuation  constant  can  be  deduced  from  the 
measured  data.  Many  optical  films  are  polycrystalline  in  nature.  The 
grain  boundaries  are  distributed  more-or-less  uniformly  throughout  the 
film,  and  act  as  the  scattering  centers.  Streaks  of  light  appear  when  light 
is  scattered  by  the  scattering  centers.  The  intensity  of  the  scattered 
light  is  proportional  to  the  intensity  of  the  guided  beam  and  can  be  used 
as  a  measure  of  the  intensity  of  the  guided  beam.  By  plotting  the  intensity 
of  scattered  light  as  a  function  of  position  along  the  propagating  beam,  one 
can  estimate  the  attenuation  constant.  A  simple  way  to  perform  the  experi¬ 
ment  is  to  use  a  TV  camera.  The  video  image  recorded  by  the  camera  can  be 
processed  electronically.  A  microscope  objective  or  an  optical  fiber  may 
be  used  to  pick  up  the  scattered  light  [13,  14].  Because  of  its  size  and 
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flexibility,  an  optical  fiber  is  preferred.  A  fixed  orientation  of  the 
fiber  relative  to  the  waveguide  and  a  constant  separation  between  the  fiber 
tip  and  the  waveguide  should  be  maintained  throughout  the  measurement. 

To  avoid  fluctuation  or  "scattering"  of  the  data,  the  numerical  aperture 
of  the  fiber  should  not  be  too  small.  The  reading  taken  when  the  fiber  is 
near  the  edges  of  the  waveguide  should  be  discarded.  It  is  also  desirable 
to  probe  the  scattered  light  as  a  function  of  position  transverse  to  the 
propagating  beam  and  to  integrate  the  signal  so  measured.  Of  course,  such 
a  measurement  would  be  very  time-consuming. 

A  Gamma  Scientific  Radiometer  interfaced  with  a  Hewlett-Packard  9825 
Computer/Controller  and  various  attachments  are  available  for  use  at  NAC. 

All  that  is  needed  is  to  construct  a  stable  platform  so  that  the  fiber  can 
be  moved  relative  to  the  waveguide  in  a  fixed,  predetermined  fashion. 

6.0  EXPERIMENTAL  OBSERVATIONS  WITH  INJECTION  LASERS 

Since  semiconductor  injection  lasers  are  small  and  can  be  modulated 
directly  at  high  frequency,  they  are  favored  in  many  communication  and 
signal  processing  applications.  The  bandgap  of  GaAs ,  or  ternary  and  quaternary 
compounds  based  on  GaAs  material,  is  such  that  radiation  emitted  from  GaAs 
injection  lasers  are  in  the  near  infrared  region  and  are  invisible  to  human 
eyes.  In  addition,  the  radiation  from  these  lasers  diverges  quickly. 

Unless  a  simple  means  is  found  to  trace  or  to  "see"  the  laser  beam,  it 
would  be  very  difficult  to  align  various  optical  components  with  respect 
to  the  invisible  beam.  Ways  must  also  be  found  to  collect  and  to  collimate 
the  beam  before  it  diverges.  Other  than  these  two  complications,  experi¬ 
ments  involving  GaAs  or  other  semiconductor  lasers  with  near  infrared  ( I R) 


emission  are  essentially  the  same  as  the  experiments  using  lasers  in  the 
visible  spectra. 

Two  simpleways  of  "seeing"  emission  in  near  IR  spectra  are  IR- 
phosphor  plates  and  image  converters.  While  the  phosphor  plates  are  easy  to 
use,  their  sensitivity  is  rather  low.  Two  crucial  components  of  an  image 
converter  are  a  photocathode  and  a  fluorescent  screen.  Electrons  are  ejected 
from  the  photocathode  as  a  result  of  the  photon  bombardment.  The  electrons 
so  released  are  accelerated  toward  the  fluorescent  screen  which  emits  light 
in  the  visible  region.  Thus  the  IR  or  near  IR  image  is  converted  into  a 
visible  one.  Depending  on  the  photocathodes  used,  the  image  converter 
may  be  used  in  different  spectral  regions.  For  example,  S-l  cathodes,  which 
are  typically  used  in  many  commercially  available  image  converters  or  IR 
viewers,  are  useful  in  the  range  between  0.3  to  1.2  pm  with  its  peak  response 
in  the  0.7  to  0.9  pm  range,  which  coincides  with  the  emission  spectra  of 
GaAs  lasers. 

TV  cameras  equipped  with  Si  cathodes  can  also  be  used  to  view  near 
IR  or  IR  emission.  However,  TV  cameras  are  too  heavy  and  bulky  for  initial 
alignment  purpose,  although  they  are  very  useful  for  the  final  viewing 
or  detailed  measurement. 

The  active,  light-emitting  area  of  an  injection  laser  is  quite  small 
and  generally  has  a  rectangular  shape.  As  a  result,  the  beam  emitted  by 
an  injection  laser  has  an  elliptical  cross  section  with  a  large  angle  of 
divergence.  Typically,  the  beam  divergences  are  10. °  (parallel  to  the  pn 
junction)  x  30.°  (perpendicular  to  the  junction),  i.e.,  170  x  520  mrad. 

For  comparison,  it  is  noted  that  the  beam  divergence  of  a  typical  HeNe 


laser  is  about  1  to  2  mrad.  Since  the  beam  diverges  fast,  a  lens  with  a 
small  F-number  or  a  microscope  objective  with  a  large  numerical  aperture 
(N.A.)  can  be  used  to  collect  the  beam  before  it  diverges.  The  F-number  of 
a  lens  is  defined  as  the  ratio  of  the  effective  focai  length  f  to  the 
diameter  d  of  the  aperture: 

F  -  number  =  f/d  ( 1 8) 


For  a  microscope  objective  situated  in  the  air  and  accepting  beams  within 
a  cone  of  half-angle  e  ,  the  numerical  aperture  is 


N.A.  *  sin  e 


(19) 


For  lenses  corrected  for  coma  and  spherical  aberration  and  with  infinite 
object  distance,  these  two  quantities  are  related: 


F  -  number  s 


1 

2N.A. 


(20) 


For  reference  purpose  and  for  convenience  of  alignment,  it  is  convenient 
to  superimpose  a  visible  light  with  the  invisible  beam.  A  mirror  and  a 
cube  beam  splitter  have  been  used  for  this  purpose. 


7.0  FABRICATION  OF  OPTICAL  WAVEGUIDES 


Many  schemes  are  available  for  fabricating  the  waveguides.  Depending 
on  the  fabrication  process,  the  index  profile  in  the  substrate  and  film 
regions  may  vary  considerably.  A  three-layer  structure,  with  nf  >  ns  >  nc» 
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?s  the  simplest  optical  waveguide  configuration.  In  most  optical  waveguides, 
the  transition  from  the  film  region  to  the  substrate  region  is  gradual  rather 
than  abrupt.  Extensive  reviews  of  all  existing  fabrication  processes  have 
been  reported  by  Chang  et.  al.t  [16],  by  Tien  et.  a!.,  [17,  18],  and  by 
Zernike  [19].  Not  all  waveguides  are  applicable  in  practical  systems,  however. 
Here,  the  fabrication  of  waveguides  which  are  potentially  useful  in  practical 
systems  are  reviewed.  Specifically,  techniques  useful  in  forming  solution- 
deposited  waveguides,  and  waveguides  based  on  glass  and  lithium  niobate 
substrates  are  considered.  In  terms  of  practicality,  waveguides  based  on 
GaAs  should  also  be  included  in  the  list.  As  is  well  known,  fabrication  of 
semiconductor  thin-films  is  a  precise  and  complex  discipline.  Numerous 
articles  and  books  are  devoted  to  the  subject  of  semiconductor  thin-films 
in  general,  and  GaAs  or  ZnO  material  [20-23]  in  particular.  Therefore, 
optical  films  based  on  GaAs  or  ZnO  materials  are  not  included  in  the  review. 

7.1  Solution  deposited  waveguides: 

Of  all  techniques  available  for  producing  optical  films  and  components, 
techniques  based  on  solution  deposition  are  the  simplest.  Most  thin-film 
fabrication  processes  have  their  origin  in  microelectronics.  These  techniques 
are  eminently  suited  for  fabricating  small  or  miniature  components.  If 
optical  films  or  components  with  large  areas  are  needed,  solution  deposition 
or  rf  discharge  polymerization  are  probably  the  best  candidates  [24-30]. 

The  properties  of  polymeric  thin-films  have  been  studied  by  Swalen  et.  al., 
[24].  Some  details  of  fabricating  photoresists  (Kodak’s  KPR,  and  Shiple 
AZ1350),  polyurethane ,  epoxy,  lead  silica,  polystyrene  (PS) ,  polymethyl 
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methacrylate  (PMMA),  styrene  acrylonitrile  copolymer  (SAN)  or  combinations 
of  PMMA  and  SAN  have  been  reported  [25  -  293 •  The  combination  of  PHMA 
and  SAN  has  an  interesting  property  in  that  the  index  of  the  film  can  be  adjust¬ 
ed  from  1 . 489  to  1.5&3  by  properly  mixing  PMMA  with  SAN.  The  pertinent  fabrica¬ 
tion  parameters  are  summarized  in  Table  1.  The  shortcomings  of  these  films 
include  a  short  life  time  and  in  some  cases,  attenuation  is  sensitive  to 
ambient  humidity  [29]. 

After  the  solution  is  prepared,  it  can  be  applied  onto  a  solid  substrate, 
usually  a  glass  slide.  Uniform  coating  can  be  achieved  by  dipping,  horizontal 
flow,  spinning  coating  or  "doctor  blading"  methods  [2^].  "Dipping"  simply 
means  dipping  the  glass  slide  into  the  solution  and  then  withdrawing  the 
slide  from  the  solution  with  constant  speed.  The  film  thickness  is  controll¬ 
ed  by  the  viscosity  of  the  solution,  and,  to  some  extent,  the  speed  of  with¬ 
drawing.  The  slide  is  dried  and/or  baked  in  a  horizontal  position.  The 
horizontal  flow  method  starts  with  covering  the  slide  with  the  solution 
from  a  syringe, the  slide  is  then  brought  to  a  vertical  position  so  that  the 
excess  solution  is  drained.  The  slide  is  returned  to  the  horizontal  posi¬ 
tion  for  drying  and/or  baking.  In  spin  coating,  the  slide  is  placed  in  a 
spinner.  After  the  slide  is  covered  with  the  solution,  it  is  spun  with  a 
specific  speed  and  for  a  specific  duration  so  as  to  achieve  a  given  thickness. 
Since  the  speed  and  duration  of  spin  can  be  controlled,  in  addition  to  the 
density  and  the  viscosity  of  the  solution,  the  films  so  obtained  are  quite 
uniform.  In  the  "doctor  blading"  method,  a  knife  edge  is  used  to  spread 
the  solution  along  the  surface  of  the  slide. 


7.2  Glass  waveguides: 

Host  optical  systems  used  in  communication  or  sensor  applications  will 
require  optical  fibers  as  their  transmission  media.  At  present,  optical 
fibers  are  made  of  glass  or  quartz  and  have  indices  of  refraction  in  the 
range  of  1.5  to  1.6.  Since  the  10  components  are  a  part  of  the  communica¬ 
tion  or  sensor  system,  it  Is  necessary  to  interface  the  10  components  to  the 
optical  fibers.  To  reduce  the  reflection  loss  at  the  interfaces,  the  indices 
of  the  glass  fibers  and  the  10  components  should  be  matched  as  closely  as 
practical.  Waveguides  based  on  glass  or  quartz  substrate  are  the  natural 
choice  in  those  applications.  R.F.  sputtering  [31  -  38],  chemical 
vapor  deposition  (CVD)  [39],  vacuum  or  electron  beam  deposition  [35,  40]  , 
ion  implantation  [41  -  43],  and  ion  exchange  [44  -  52]  methods  have  been 
used  in  forming  waveguides  on  various  glass  slides.  Thermal  evaporation 
in  vacuum  is  probably  the  simplest  way  of  depositing  thin  films  onto  a 

substrate.  However,  films  produced  by  thermal  evaporation  tend  to  be 

\ 

quite  lossy.  If  the  electron-beam  evaporation  process  is  used,  films 
2  to  A  pm  thick  may  be  obtained  with  an  attenuation  slightly  larger  than  1.0  dB/cm 
[40]  which  is  comparable  to  the  waveguides  made  by  sputtering  processes. 

Films  of  good  quality  with  various  compositions  may  be  produced  by  the 
r.f.  sputtering  process.  The  composition  of  the  gases  in  which  sputtering 
is  performed,  the  rf  power  and  duration  of  sputtering  can  be  controlled  with 
ease.  The  gas  pressure  and  the  temperature  of  the  substrate  can  also  be 
adjusted  to  some  extent.  These  fabrication  parameters  can  be  varied  to 
form  films  having  the  desired  composition  and  thickness.  Obviously,  the 
film  thickness  is  directly  proportional  to  the  duration  of  sputtering.  The 
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rf  power  can  be  varied  to  tailor  the  index  of  sputtered  film,  and  the 
stoichiometry  of  the  film  can  be  compensated  by  admitting  various  gases  in 
the  various  stages  of  sputtering  [33,  3*0.  The  sputtering  parameters  and 
the  properties  of  the  films  so  obtained  are  summarized  in  Table  2. 

The  refractive  index  of  a  film  is  a  function  of  the  density  and  the 
electronic  polari zabi 1 i ty  of  the  constituent  species,  and  is  described  by 
the  Claus i us-Hossotti  relation  [53].  If  the  density  of  the  material  is 
increased,  so  is  the  index  of  refraction.  More  importantly,  if  the  composi¬ 
tion  of  the  material  is  changed  so  that  species  with  small  electronic 
polarizability  are  replaced  by  species  with  large  electronic  polarizability, 
the  index  of  refraction  is  also  increased.  In  forming  waveguides  by  solution 
deposition,  thermal  evaporation  or  r.f.  sputtering,  layers  of  material  with 
higher  index  are  depos i ted  onto  substrates,  and  the  resulting  index  profiles 
are  more-or-less  abrupt.  In  the  ion  exchange  process,  however,  an  optical 
waveguide  is  formed  when  ions  with  small  pol a r i zab i 1 i ty  near  the  substrate 
surface  are  replaced  by  ions  with  larger  polarizability,  and  the  index 
profile  of  an  ion-exchanged  waveguide  changes  gradually  from  the  surface 
(i.e.,  the  "f i lmu  region)  to  the  substrate  region.  Ion  exchange  processes 
can  be  used  to  form  waveguides  in  glass,  lithium  niobate  or  other  substrates. 
In  this  section,  the  discussion  is  restricted  to  glass  or  quartz.  Typically, 
a  glass  slide  is  immersed  at  elevated  temperature  in  molten  salts  containing 
the  ions  to  be  exchanged.  Compounds  such  as  KNO^  ,  TINO^,  and  NaNO^  [44],  or 
AgNOj  [45”47,  50], or  AgNO^  and  NaNO^  [51],  or  a  eutectic  melt  of  L^SO^ 
and  I^SO^  [48]  have  been  used  in  conjunction  with  various  glasses. 
Alternatively,  a  solid  silver  film  or  a  vapor  stream  of  silver,  produced 


by  thermal  evaporation,  may  be  used  instead  of  melts  [49,  52].  The  exchange 
process  can  be  done  with  or  without  the  assistance  of  an  electric  field. 

In  the  presence  of  an  electric  field,  the  speed  of  ion  exchange  is  increased, 
and  processing  time  can  be  reduced.  Electric  fields  can  also  be  used  to 
make  the  index  profile  more  abrupt  and  to  move  the  peak  of  the  index  profile 
to  the  interior  region  of  the  substrate.  By  diffusing  Ag+,  Tl+  and  K+  ions 
into  the  glass  to  replace  Na+  and  Li+  ions,  the  index  of  the  diffused  region 
will  be  higher  than  the  index  of  the  undiffused  region.  Masks  with  windows 
or  apertures  can  be  placed  on  the  glass  surface  so  that  ion  exchange  is 
restricted  to  the  selected  regions.  In  this  fashion,  strip  waveguides  can 
be  formed  [46].  The  fabrication  parameters  of  ion  exchange  process  are 
tabulated  in  Table  3- 

7.3  Lithium  niobate  waveguides: 

Since  lithium  niobate  and  lithium  tantalate  are  highly  transparent 
and  have  large  electro-optic  coefficients,  they  are  commonly  used  waveguide 
materials.  In  particular,  lithium  niobate  waveguides  have  been  used  to 
demonstrate  the  feasibility  of  many  device  concepts  and  practical  systems. 

The  epitaxial  growth  by  melting  method  has  been  used  to  grow  a  solid  solu¬ 
tion  of  (LiNbO^)a(LiTaOj) |_a  onto  LiTaOj  substrate  [54].  A  guiding  layer 
can  also  be  formed  by  thermally  diffusing,  or  diffusing  under  electric 
field,  metals  like  Cu,  Al ,  Ge,  Cr,  Fe,  and  Nb  into  LiTaO^  [55  "  57],  or 
by  releasing  L^O  from  the  surface  of  LiNbOj  or  LiTaOj  [58,591*  By  far 
the  most  popular  method  is  to  deposit  a  layer  of  transition  metal  like  Ti , 

V  or  Ni,  on  the  surface  of  LiNbO^  and  to  diffuse  these  elements  into  the 


substrate.  This  is  known  as  the  in-diffusion  process  [60] .  Of  the  transi¬ 
tion  metals  mentioned,  Ti  is  the  most  popular  element  used,  and  it  leads 
to  the  well  known  T i : L t NbO^  waveguides  [60  -  64].  In  the  diffusion  process, 
the  temperature  of  the  substrate  has  to  be  raised,  l^O  is  released  from  the 
surface  resulting  in  an  unwanted  out-diffused  waveguide  at  the  surface 
of  the  substrate.  Considerable  effort  was  devoted  to  the  elimination  of 
the  out-diffused  layer  [65  -  70] .  The  simplest  way  of  suppressing  the  out- 
diffused  waveguide  layer  is  to  add  water  vapor  in  the  diffusion  gas  [70]. 
Table  4  lists  most  pertinent  parameters  in  forming  in-diffused  waveguides. 

Ion  exchange  methods  have  also  been  used  to  form  optical  waveguides 
on  LiNbOj  and  LiTaO^  [71,  72].  By  exchanging  Li+  ions  or  the  vacant  sites 
for  Li  +  ions  with  metallic  ions,  optical  waveguides  can  also  be  formed 
in  LiNbOj  of  LiTaO^.  By  placing  polished  x-cut  plate  of  LiNbO^  in  molten 
AgNO^  at  360°C  for  3  to  15  hours,  waveguides  capable  of  supporting  1  to 
3  TE  modes  were  obtained  by  Shah  [71].  The  attenuation  constant  was  measured 
with  an  optical  fiber  probe  to  be  6  dB/cm.  The  increase  in  the  extraordinary 
index  of  refraction  is  quite  large,  (Ane  z  0.12).  But  the  change  in  the 
ordinary  index  of  refraction  is  insignificant,  and  no  TM  modes  were  observed. 
Experiments  were  also  performed  showing  that  the  electro-optic  coefficients 
of  LiNbOj  are  not  significantly  affected  by  the  ion  exchange  process.  Y-cut 
LiNbO^  or  LiTaO^  samples  were  also  placed  in  AgNO^  melt  for  24  hours,  no 
waveguide  layer  was  formed,  however.  Additional  work  reported  by  Jackel  [72] 
showed  that  the  loss  can  be  reduced  to  the  0.8  to  2.0  dB/cm  region.  Larger 
change  in  the  extraordinary  index  of  refraction  has  been  obtained  when  the 
Li+  ions  are  exchanged  with  Tl+  ions  in  molten  TINO^. 


8.0  CONCLUSION 


With  the  introduction  of  integrated  optics,  a  new  class  of  optical 
components  and  systems  emerge.  Thin-film  optical  waveguides  are  the  building 
blocks  of  these  systems.  Methods  for  fabricating,  and  experimental  and 
numerical  techniques  for  characterizing  these  waveguides  were  studied. 
Techniques  useful  for  performing  experiments  in  visible  and  infrared 
spectra  were  discussed.  This  study  serves  as  a  brief  introduction  to  a 
new  horizon  in  command,  control  and  communication  systems  which  has  particular 
utility  in  avionics  instrumentation. 
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APPENDIX  -  LISTING  OF  A  COMPUTER  PROGRAM 


In  this  Appendix,  the  program  MAIN  and  the  subroutines  used  to  compute 
or  estimate  the  refractive  index  n^  and  the  thickness  h  of  the  film  from  the 
measured  coupling  angles  are  listed.  These  programs,  written  in  the  standard 
Fortran  language,  have  been  tested  as  described  in  the  main  text  of  this 
report . 


program  main ( input;  output ) 

common  /dl/  lambda#  nefms (20) #  m<20)#  nss#  ncs#  mrho#  nlines 

common  /d3/  neft(20) 

real  ns#  nc#  nf  #  np#  lambda 

real  nss#  ncs#  nfs#  kw 

dimension  theta(20)#  m<20) 

real  neft(20)#  nefm<20),  nefms(20) 

external  ersum 

open ( 1  #  file='mdata'#  status='old') 
rewind  1 

r ead  ( 5#  180 ) ns #  nc . np #  lambda#  ep sm#  mr h o 

format  ( 5<  f  8-  4#2x),  i2) 

anp2=np*np 

nss=ns*ns 

ncs=nc*nc 

p i=3.  141592654 

if(mrho  .  eq.  0)  write  (6# 190) 

if(mrho  .  eq.  1)  write  (6#  191) 

format <3x,  'TE  MODES'  ) 

format (3x, 'TM  MODES'  ) 

write (6#  183 ) n s #  nc #  np #  ep sm,  lambda 

forma t  < 2x *  'ns-  'f8.  4#  '  nc=  'f8.  4#  '  np=  'f8.  4#  9  angle  oi 

1  f 9.  4#  3x  i  '  lambda*  'f8.  4) 
epsmr  =  ep  sm*p i /180- 
do  184  i  =  1 #  20 

read ( 5# *#  end =299)  m(i)*theta(i) 
nl ines=i 


np=  'f8.  4#  9  angle  of  prism  * 


write (6#  185) 

format(3x#  'More  than  20  m-lines  ? 
goto  999 

do  298  i=l, nlines 

write(6# 186)  m<i>#  theta(i) 

f  ormat  <3x,  i2#  f  8.  4,  '  degrees') 


c 
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do  220  ie*l.nlines 
thetar*theta(  ie)*pi/180. 

nefm( ie )*s in< the tar )#co» ( epsmr >+( tin <  epsmr )*sqrt ( anp2— sin ( thetar ) 
1  *sin(thetar) ) ) 
nefms  (  ie)*nefm(  ie )*nefm( ie  > 
nef t ( ie )*nefm( ie ) 

write (6. 215)  ie. theta ( ie )# nefm( ie) 

215  formatOx.  i2.  '-th  coupling  angle  *  'f8.  3*3x«  'N  eff  ■  'f8.  4) 

220  continue 

c  if(nlines  .  gt.  2)  goto  500 

nf*l.  048*amaxl (nefm( 1 )»  nefm(2)) 
nf s*nf *nf 
fts-1. 
f  tc*l. 

do  300  iq-1.  15 

if(mrho  .  eq.  1)  f tc*(nf s/ncs >*(nf s/ncs) 
if(mrho  . eq.  1)  f ts=(nfs/nss)*(nf s/nss) 
phis*atan(sqrt(fts«(nefms( 1 )  -  nss)/(nfs  -  nefms(l)))) 
phic=atan(sqrt(ftc*(nefms(l)  -  ncs)/(nfs  -  nefms(l)))) 
psi l=m( 1 )*p i+ph is+ph ic 

ph is-atan ( sqrt ( f ts*(nefms <2)  -  nss)/(nfs  -  nefms(2)))) 
ph ic-atan ( sqrt ( f tc*(nefms (2)  -  ncs)/(nf*  -  nefms(2)))) 
psi2=m(2)*p  i+ph  is+p.hic 

f* (nefms (2)*psi l*psi 1 -nefms ( 1 )*psi2*psi2)/ (psil*psil  -psi2*psi2) 
test=abs( (nfs-f )/nfs) 
ff*sqrt(f > 

write(6,  315)ff»nf»  test 

315  formatOx.  'nf  =  'el4.  6»3x.  'last  nf=  'e!4.  6/  'test=  '  el4.  6) 
if  (test  .It.  0.000001)  goto  400 
nf*ff 
nf  s=f 

300  continue 

write (6.  320) 

320  formatOx.  '***  2  m-lines.  not  converging  in  15  iterations.***') 
goto  999 

400  wl*lambda*psi  1/  (2.  *p i*sqr t ( nf s-nef ms ( 1 ) ) ) 
w2*lambda*psi2/ (2.  *p i*sqrt (nf s-nefms (2) ) ) 
write(6, 316)nf.  wl.  w2 

31£  formatOx.  'nf*  'f8.  4.  3x,  'wl*  'f  10.  4,  3x.  'w2*  'fl0.4,3x,  '  microns') 
if(nlines  . eq.  2)  goto  999 
500  kw*2.  *p i*wl/lambda 

hx=l.  e-8 
hy*i.  e  -8 
epsiln*l.  e-6 

call  grdt( er sum. hx>  hy .  nf.  kw.  epsiln) 
w=kw*lambda/ (2.  *pi) 
write(6. 349>nf»  w 

349  format(3x>  'nf*  'f8.  4.3x.  'w=  'flO.  4.3x.  '  microns') 

999  stop 


subroutine  disp(aneffs#  dispft) 

common  /dl/  lambda#  nefms(20)#  m(20)*nss#ncs#mrho#nlines 

common  /d2/  nfx#  uig#  mi  j#  factor 

real  nf#nfx#nss#ncs#nfs 

dimension  m(20) 

nf s*nf  x*nf  x 

pi =3.  141592654 

anef s*anef f s/fac  tor 

fts-1. 

ftc  =  l. 

if(mrho  .  eq.  1)  f tc*(nf s/ncs )#(nfs/ncs) 
if(mrho  .  eq..  1)  f  ts® <nf  s/nss  )*(nf  s/nss ) 
ph is=atan ( sqrt  <  f ts*(anef s-nss > / (nf s  -anef s ) ) > 
ph ic®a tan ( sqrt ( f tc*(anef s-ncs ) / (nf s  -anefs) ) ) 
psi=mi  j*p  i+ph  is+ph  ic 
d  ispf  t=wy*sqrt  (nf  s-anef  s  >-  psi 
write (6. 980 )nf x» wg#aneffs#  dispft 
980  format  (lx#  'disp'#4(e!4.  6#lx>) 
return 
end 


function  ersum(rtn. rtw) 
common  /dl/  lambda* nefms (20) «  m(20> * nss. ncs* mrho. nl ines 
common  /d2/  nfx*  wy.  mi  j*  factor 
common  /d3/  ncft(20) 
real  nef t (20) .  nefmi j 
real  lambda*  nefms (20) * nss. ncs. nf x 
dimension  m(20> 
complex  rtsf(4) 
external  disp 
logical  fnreal 
real  nfse 
kn*0 
n=l 

epl*l.  e-6 
ep2»l.  e-6 
fnreal*.  true, 
max  its 100 
ersum*0. 
factor=l.  e+4 
nf x=rtn 
Uiyartlll 
nf se=nf x*nf x 
i jn=nss*f actor 
i jm=nfse*f actor  -1 
do  100  ij=l*nlines 
i jkm=0 
mi j=m( i j ) 

nefmi j*sqrt (nefms ( i j ) ) 
rts*(neft(i j) >*(neft(ij) ) 
rtsf ( 1 )*rts*f actor 

call  muller(kn*  n>  rtsf*  ma  x  i  t.  epl. ep2. disp*  fnreal ) 
call  d  isp (rtsf ( 1 > *  dispft) 
r ts*rea l(rtsf(l)/factor) 

if(rts  .  gt.  nss  .and.  rts  .It.  nfse)  goto  220 
if(ijkm  .gt.  0)  goto  250 
i jkm=l 

ansf *nss*f actor 
call  d i sp (ansf * f tr 1 ) 
do  240  ijk=ijn»  ijm*2 


do  240  ijk=ijn»  ijm*2 
ansf=i jk 

call  d isp (ansf • ftr ) 
if (ftr*ftrl  .  ge.  0.  )  goto  249 
r tsf ( 1 )=  ansf 
goto  210 
f tr l=f tr 

wr ite (6> 251 ) i jk*  rtn* rtw.  ansf*  ftr 
format  ( lx*  'ersum  ' i3*  2(  1  x«  f  10.  7).  lx.  'ansfj 
1/  4x.  'but.  no  root  is  found.  '  ) 
continue 

neft(ij)*sqrt(rts> 
i jkm*0 

ersum=ersum  +  (rts  -nefmi j )*(rts-nefmij ) 

write  (6.  998)  i  j.  rtn.  rtw.  rts.  dispft 

format  (lx.  'ersum  'i3.2(lx.fl0.  7).  lx.  'rts* 

continue 

return 

end 


'el4.  6. 


'el 4.  6.  lx.  'disp*  'el 4.  6) 
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subroutine  grdt(fn*hx»hy«rtx,rtyiep»iln) 

external  fn 

noit»0 

100  noit=noit+l 

fO*=Fn(rtxi  rty ) 
f l*f n ( rtx+h x<  rty) 
f2**fn(rtx»  rty+hy ) 
n (rtx-h x.  rty) 

•f4«fn(rtx  /  rty-hy ) 
f5*fn(rtx+hx«  rty+hy) 
f6*fn(rtx-hx»  rty+hy) 
f7=fn(rtx-hx.  rty-hy) 
fS^fnCrtx+hx.  rty-hy) 
sx=(f l-f3>/(2.  *h x ) 
sy  =  (f 2-f 4) / (2.  *hy ) 

*xx=(fl  -2.  *f 0  +f3)/(hx*hx) 
syy*=(f2  -2.  0  +  f4)/(hy*hy) 

sxy  =  ( f 5  -  -f 8  +  f 7) / (4.  *h x*hy ) 

write  (6.  141  )fO.  f  1.  f2»  f3.  f4«  F5»  #6.  f7,  F8«  sx.sy.sxx.syy.exy 

141  format(3x»  'grdt  fO»  f  1.  f2#  f3.  f4  f5.  f6.  F7«  f8  sx.  sy  sxx»syy>sxy  ' 
1/3 x,  e  12.  4.  /4(3x»  el2.  4)/4<3x,  el2.  4)/4(3x»  el2.  4)/2(3x.  el2.  4) 
2/3(3x >  el2.  4) ) 

det=s x x*sy y  -sxy*sxy 
rtxl=rtx  +  (sy*sxy  -sx*syy)/det 
rtyl=rty  +  <sx*syy  -sy**xx)/det 
write(6»  142)det#  rtx«  rty /  rtx  1*  rty  1 

142  FormatOx,  'det*  'el4.  6> 

l/3x.  'rts<  rty=  '2(3x.  el4.  6)  /3x,  'rt x 1 , rty 1  =  '2(3x,  el4.  6) ) 
error=sqrt( (rtxl-rtx)*(rtxl-rtx)  +(rtyl-rty)*(rtyl-rty)) 
error=error/sqrt<rtx*rtx  +rty*rty) 
write (6/ 998) no it 
if(error  . le.  epsiln)  goto  110 
r t  x=r  t  x 1 
rty=r ty 1 

998  •FormatOx/  'grdt  noit*  'i4> 
noit=noit+l 

if(noit  .  gt.  20)  write(6<  120) 
if(noit  .  le.  20)  goto  100 

120  format<3x»  'more  than  20  iteration.  ') 

110  return 
end 
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Figure  6  Transformation  of  Gaussian  beam  by  a  lens 
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Figure  7  Top  vieu  of  a  parabolic  waveguide  horn  (after  Burns/  et.  al.) 
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